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ABSTRACT 

The Carina dwarf spheroidal galaxy is the only galaxy of this type that shows clearly episodic star 
formation separated by long pauses. Here we present metallicities for 437 radial velocity members of this 
Galactic satellite. The metallicities and radial velocities were measured as part of a Large Programme 
with the Very Large Telescope at the European Southern Observatory, Chile. We obtained medium- 
resolution spectroscopy with the multi-object spectrograph FLAMES. Our target red giants cover the 
entire projected surface area of Carina. Our spectra are centered at the near-infrared Can triplet, which 
is a well-established metallicity indicator for old and intermediate-age red giants. The resulting data 
sample provides the largest collection of spectroscopically derived metallicities for a Local Group dwarf 
spheroidal to date. Four of our likely radial velocity members of Carina lie outside of this galaxy's 
nominal tidal radius, supporting earlier claims of the possible existence of such stars beyond the main 
body of Carina. We find a mean metallicity of [Fe/H] ~ — 1.7dex on the metallicity scale of Carretta 
& Gratton (1997) for Carina. The formal full width at half maximum of the metallicity distribution 
function is 0.92 dex, while the full range of metallicities is found to span ~ —3.0 < [Fe/H] < 0.0 dex. 
The metallicity distribution function might be indicative of several subpopulations distinct in metallicity. 
There appears to be a mild radial gradient such that more metal-rich populations are more centrally 
concentrated, matching a similar trend for an increasing fraction of intermediate-age stars (Harbeck et 
al. 2001). This as well as the photometric colors of the more metal-rich red giants suggest that Carina 
exhibits an age-metallicity relation. Indeed the age-metallicity degeneracy seems to conspire to form 
a narrow red giant branch despite the considerable spread in metallicity and wide range of ages. The 
metallicity distribution function is not well-matched by a simple closed-box model of chemical evolution. 
Qualitatively better matches are obtained by chemical models that take into account also infall and 
outflows. A G-dwarf problem remains for all of these models. 

Subject headings: Galaxies: abundances — Galaxies: dwarf — Galaxies: evolution — Galaxies: stellar 
content — Galaxies: structure — Galaxies: individual (Carina) — Local Group 



1. INTRODUCTION 

Dwarf galaxies come in many flavors ranging from gas- 
rich to gas-poor, from irregular in shape to ellipsoidal, 
from actively star-forming to quiescent (Grebel 2001). The 
least massive, least luminous galaxies known are the dwarf 
spheroidal (dSph) galaxies. They are characterized by 
absolute V-band luminosities My > —14 mag, surface 
brightnesses of /xy > 22 mag arcsec"^, Hi masses of 
A^Hi ^ IO^Mq, estimated total masses Mtot of a few times 
10 ''Mq, have shallow projected radial light profiles and 
tend not to be rotationally supported (see Grebel, Gal- 
lagher, & Harbeck 2003 and references therein). Because 
of their high velocity dispersion, their velocity dispersion 
profiles, their morphology and their lack of depth extent 
many dSphs are believed to be dominated by dark mat- 
ter (e.g., Mateo 1997; Odenkirchen et al. 2001; Klessen, 
Grebel, & Harbeck 2003; Wilkinson et al. 2004). In galaxy 
groups dSphs are usually found within ~ 300 kpc around 
more massive galaxies (Fig. 1 in Grebel 2005). The gas 



deficiency of dSphs remains an unsolved puzzle - dSphs 
typically contain even less gas than expected from red gi- 
ant mass loss over time scales of several Gyr. Searches 
for neutral and ionized gas usually only lead to low upper 
limits on the gas content (e.g., Gallagher et al. 2003 and 
references therein). 

DSphs are typically dominated by either very old pop- 
ulations (> 10 Gyr) or by intermediate-age populations 
(1 — 10 Gyr). Very few show more recent star forma- 
tion. All dSphs (and all other dwarfs) studied in sufficient 
detail so far have been shown to contain ancient popu- 
lations that are indistinguishable in age with the oldest 
age-datable populations in the Milky Way (Grebel 2000, 
Grebel & Gallagher 2004). The detailed star formation 
histories and the metal enrichment of dwarf galaxies in- 
cluding dSphs vary widely; no two dwarfs are alike (e.g., 
Grebel 1997, Mateo 1998). DSphs usually show contin- 
uous star formation with some amplitude variations and 
declining rates at more recent times. A number of dSphs 
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are entirely dominated by old populations, for instance, 
Draco and Ursa Minor. Younger and/or more metal-rich 
populations in dSphs are more centrally concentrated, in- 
dicating extended star formation episodes in the centers of 
their shallow potential wells (Harbcck ct al. 2001). 

The Galactic dSph Carina stands out among the dSphs 
in the Local Group because of its unusual, episodic star 
formation history. In no other dSph clear evidence for 
well-separated episodes of star formation has been found. 
Carina was discovered in 1977 (Cannon. Hawardcn. & 
Tritton 1977) on an ESO/SRC(J) Southern Sky Survey 
plate. Carina is located at a heliocentric distance of 
94 ± 5 kpc and belongs to the fainter dSphs (My = —9.4, 
Hv = 25.5 ± 0.4; see Grebel 2000 for references). First in- 
dications that Carina is not a purely old, globular-cluster- 
like system came from the discovery of carbon stars (Can- 
non, Niss, & Norgaard-Nielsen 1981; Mould et al. 1982). 
Subsequent deep photometry revealed a younger main se- 
quence (Mould & Aaronson 1983; MighcU 1990a,b) and 
showed that the bulk of the stars in Carina is of interme- 
diate age (i.e., younger than 10 Gyr). However, an old 
population (traced by a horizontal branch and RR Lyrae 
stars, e.g., Saha, Monet, & Seitzer 1986) is present as well. 
Mighell's (1990a, b) data show two distinct main-sequence 
turnoffs. Smecker-Hane et al.'s (1994) color-magnitude 
diagram (CMD) of Carina reveals a morphologically dis- 
tinct, prominent red clump of an intermediate-age popula- 
tion next to a woU-dcfincd red and blue horizontal branch. 
Mighell and Smecker-Hane et al. concluded that Carina 
must have undergone at least two separate episodes of star 
formation. 

The episodic nature of the star formation in Carina was 
then confirmed by deeper photometry: Smecker-Hane et 
al. (1996) presented a color-magnitude diagram (CMD) 
that reached i? ~ 25 mag. Three well-defined, distinct 
main-scqucncc turn-offs (MSTOs) corresponding to ages of 
2, 3-6, and 11-13 Gyr are visible in these data. Smecker- 
Hane et al. note that each of the three MSTOs con- 
nects to the same narrow red giant branch and conclude 
that "regardless of age. Carina stars are metal-poor with 
[Fc/H] ^ —1.86 dcx and a spread in metal abundance of 
o-[Fe/Hi < 0.2 dex." Mighell (1997), Hurley-Keller, Ma- 
teo, & Nemec (1998), MoneUi et al. (2003), and Rizzi 
et al. (2003) arrive at similar conclusions based on addi- 
tional deep CMDs^ Hurley-Keller et al. (1998) argue that 
only 10- 20% of Carina's stars formed very early on, while 
the majority formed approximately 7 Gyr ago and ~ 30% 
about 3 Gyr ago. (Note that differences in the derived ages 
or times of the extended episodes of star formation arc 
mainly caused by the use of different isochrone models). 
MoneUi et al. (2003) suggest that the blue plume stars of 
Carina may be as young as 1 Gyr, and that the anomalous 
Cepheids of Carina may belong to an even younger popu- 
lation with an age of ^ 0.6 Gyr (sec also Mateo, Hurley- 
Keller, & Nemec 1998; Porctti 1999; and McNamara 2000 
for the identification and the discussion of these stars). 
The younger populations of Carina are more centrally con- 
centrated than the extended distribution of the old stars 
(Harbeck et al. 2001; MoneUi et al. 2003), similar to the 
trend seen in other dSphs. 



Smecker-Hane et al. (1999) obtained Can triplet (CaT) 
metallicities for 52 red giants in Carina. These authors 
argue that the narrow red giant branch of Carina results 
from an age-metallicity conspiracy in the sense that more 
metal-rich, but yoimgcr stars come to lie at the same lo- 
cation in the color-magnitude plane as older, metal-poor 
stars, and that photometry alone underestimates the true 
metallicity spread in Carina. Based on a photometric 
study, Rizzi et al. (2003) suggest that the narrow red giant 
branch is a consequence of the contribution of the dom- 
inant intermediate- age star- formation episode, while the 
contribution of the ancient episode is almost negligible. 

All these findings underline the highly complex star for- 
mation history of Carina. There is not yet a satisfactory 
explanation why Carina would have experienced episodic 
star formation with extended quiescent episodes in be- 
tween, and why its evolution was so different from that 
of other dSphs. As pointed out by Smecker-Hane et al. 
(1996), episodic accretion of fresh gas would be difficult 
given the inferred mass of Carina's dark matter halo, and 
cooling and sinking back of previously ejected gas should 
have led to higher enrichment. 

Our current knowledge of the detailed evolutionary his- 
tory of nearby dwarf galaxies is mainly based on photom- 
etry supplemented by rather sparse spectroscopic infor- 
mation (e.g., Grebel 1999). However, this situation is 
now changing as a growing body of spectroscopic infor- 
mation is becoming available. Spectroscopy plays a par- 
ticularly important role in permitting us to break the age- 
metallicity degeneracy that plagues purely photometric 
color-magnitude diagram analyses. When the metallic- 
ity of individual stars can be measured independently via 
spectroscopy, this information can be added to constrain 
photometric determinations of the star formation history. 
For gas-deficient galaxies like dSphs, our primary source 
of metallicity information arc their red giants, which are 
now easily accessible for ground-based 8 to lOm-class tele- 
scopes. Ultimately, wc may be able to derive detailed age- 
metallicity relations. Spectroscopic information for a siz- 
able number of stars is gradually becoming available for 
an increasing number of dwarf galaxies (e.g., Cote, Oke, 
& Cohen 1999; Smecker-Hane et al. 1999; Guhathakurta, 
Reitzcl, & Grebel 2000; Tolstoy et al. 2001; Pont et al. 
2004). In the current study we add spectroscopic metallic- 
ities for several hundred red giants in the enigmatic dSph 
galaxy Carina. Individual element abundance ratios in 
stars allow one to impose strong constraints on the modes 
of star formation in galaxies. The present study only deals 
with the measurement of overall stellar metallicities. The 
discussion and analysis of elemental abundances will be 
left for a future paper (Koch et al., in prep.). 

Here we present a comprehensive spectroscopic study of 
Carina in order to analyze the overall metal content of sev- 
eral hundred red giants distributed throughout this dSph 
galaxy. In the current paper we analyze the metallicities of 
these stars in order to investigate several of the questions 
raised in the previous paragraphs. In particular, we wish 
to measure the true, spectroscopic mean metallicity and 
metallicity spread of Carina, derive its metallicity distri- 
bution function, investigate whether age and metallicity 



^ We will not discuss here star formation histories of Carina derived from HST photometry since the field coverage and the number of stars 
are too small to do so reliably. 
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do indeed conspire, search for possible spatial gradients, 
and explore its evolutionary history taking into account its 
chemical enrichment. This paper is organized as follows: 
Sect. 2 introduces our data and describes the reduction 
procedures. In Sect. 3 we present the steps taken to cali- 
brate our measured CaT equivalent widths and put them 
on a reference metallicity scale. The resulting metallicity 
distribution function (MDF) is presented in Sect. 4. In 
this Section we also discuss the implications for a possible 
age-mctallicity relation in Carina. In Sect. 5, we investi- 
gate the existence of radial gradients in the MDFs. Sect. 
6 then discusses the results in the context of simple mod- 
els of chemical evolution. In Section 7, we summarize our 
findings. 

2. OBSERVATIONS AND REDUCTION 

In the course of the ESO Large Programme 171. B- 
0520(A) (PI: Gilmore), which is dedicated to the analy- 
sis of kinematic and chemical characteristics of Galactic 
dSphs, we observed five fields in Carina. These fields were 
selected such that we cover most of the area of Carina, 
but also go a bit beyond its nominal tidal radius. In these 
fields we obtained spectra of red giants with the goal of de- 
termining their metallicities through the well-established 
CaT method and of measuring their radial velocities to 
constrain the dark matter profile of Carina. In addition 
four globular clusters (NGC3201, 4147, 4590 and 5904) 
were observed in order to permit us to place our CaT 
measurements on a scale of known reference metallicities 
(Rutledge et al. 1997a). 

2.1. Target selection 

Targets in Carina were chosen from photometry and as- 
tromctry obtained by the ESO Imaging Survey (EIS; sec 
Nonino et al. 1999 for details) during the Pre-FLAMES 
era using the Wide Field Imager (WFI) at the 2.2 m tele- 
scope at La Silla, Chile. In the framework of EIS a num- 
ber of different targets and fields was observed, includ- 
ing a region of four WFI fields, comprising approximately 
1 deg^, located at the center of the Carina dSph. The re- 
duced and transformed photometry is available from the 
EIS web pages at http://www.eso.org/scicnce/eis/. For 
our spectroscopic study, we need to pick sufficiently bright 
objects. Moreover, the CaT method is only calibrated for 
intermediate-age and old red giants (Cole et al. 2004), so 
we chose luminous red giant branch (RGB) stars for the 
present work. 

We selected our targets to cover magnitudes ranging 
from the tip of the RGB down to 3 mag below the RGB 
tip, thus down to 20.3 mag in V-band apparent magnitude. 
This corresponds to an absolute magnitude My » 0.3 at an 
adopted apparent V-band distance modulus of 20.05 and 
a reddening of 0.05 (Mighell 1997). The actual values of 
the distance and reddening are not critical for the analysis 
in this paper, since we utilize apparent magnitudes rela- 
tive to the horizontal branch (HB) locus. This differential 
approach also reduces the problems inherent in the EIS 
photometry, which go back to intrinsic problems of WFI 
photometry (see, e.g., Koch et al. 2004a,b). Our chosen 
RGB magnitude range is such that even for the faintest 
stars one is generally able to achieve high signal-to-noise 

^ For details on the individual components of the FLAMES spectrograph 



ratios (S/N ~ 25) within a reasonable integration time 
and under good observing conditions (see below). Further- 
more, we selected stars across the full width of the RGB 
(approximately 0.2 mag in B — V) in order to ensure that 
we were unbiased with respect to metallicity or age, and 
include also potential extremely metal-poor and metal-rich 
giants (see the CMD in Fig. 1). With these constraints and 
with the limitations in possible positions on the fiber plate 
in order to avoid fiber crossings, a total number of 1257 
stars was targeted. However, our photometric selection is 
inevitably contaminated by Galactic foreground stars. 

For our calibration clusters, we used the standard fields 
of Stetson (2002). More than 80 red giant candidates al- 
together were picked from Stetson's B- and V-band pho- 
tometry. These clusters range in metallicity from approx- 
imately — 1.1 dex to ~ —2 dex in [Fe/H] (see Table 1 for 
details). 

The available photometry of Carina targets is, in con- 
trast to the Stetson list, in the Johnson-Cousins system 
provided by the EIS. Because of differences in the filter 
curves, this system differs from the standard Johnson- 
Cousin system as defined by Landolt's (1992) and Stet- 
son's (2000) UBVRI standard stars. Conveniently, the 
observed EIS fields overlap with a standard star field in 
Carina (Stetson 2000) allowing a direct comparison be- 
tween the ElS-photometry and that provided by Stetson. 
We obtained a linear transformation using the 258 stars 
in common, and used this to place both the calibration 
clusters and the science targets on Stetson's homogeneous 
photometric standard system, avoiding any offsets between 
the different sets of photometry. 

2.2. Data acquisition 

Our observations were performed in visitor mode with 
the FLAMES (Fibre Large Array Multi Element Spectro- 
graph) facility at the Very Large Telescope UT2 of the Eu- 
ropean Southern Observatory on Cerro Paranal in Chile. 
We observed during 22 nights spread over two semesters 
in 2003 and 2004 (as detailed in Table 2). We used the 
GIRAFFE spectrograph with both MEDUSA^ fiber sHts 
in "low-resolution" mode. We observed with grating L8 
centered at the near-infrared CaT. This set-up provides 
a nominal resolving power of R=6500 at the center of 
the spectra and covers the wavelength range 8206-9400 A. 
Each pointing is designed to handle 132 fibers, distributed 
across a field of view of 25' diameter. We dedicated about 
20 fibers per configuration to blank sky for sky subtraction 
purposes. 

We observed five different fields in and around the cen- 
ter of Carina (Table 3), which we observed with several 
configurations of the fiber positioning plates in order to 
maximize the number of targeted stars. The location of 
our fields is displayed in Fig. 2. The total required ex- 
posure time per configuration in Carina was six hours in 
order to reach a nominal S/N of at least 20. This S/N is 
needed to obtain highly accurate equivalent width (EW) 
measurements at our spectral resolution. The observations 
were split into single exposures of approximately one hour 
each to facilitate cosmic ray removal. 

Afost of the nights were clear and good sky conditions 
resulted in a typical seeing of 0'.'5-0'.'7, although the first 

see http://www.eso.org/instruments/FLAMES/ 



4 



Koch et al. 



run in 2003 February was hampered by non-photometric 
conditions and seeing reached values as bad as 2". Spec- 
tra with low S/N ratios were excluded from the chemical 
abundance analysis, where the lower limit for marginally 
measuring EWs is at w 10, but were still useful for the 
kinematic analysis (Wilkinson et al., in prep.). 

2.3. Data reduction 

In the course of ESO's standard calibration exposures, 
bias frames, fiber flatfield exposures and Th-Ar spectra 
from a calibration lamp were taken automatically during 
daytime. We reduced our data using version 1.09 of the 
FLAMES data-reduction system, girbldrs, and the respec- 
tive pipeline girbldrs-pipe-1.05 (Blecha et al. 2000). The 
basic steps comprise bias correction as well as fiber local- 
ization and successive adjustment via the flatfield expo- 
sures. 

Our spectra were extracted by summing up the pixel 
values along a virtual slit with a width of one pixel. The 
final re-binning to linear wavelength space was obtained 
with the master wavelength solution from the calibration 
Th-Ar spectra, which were reduced and extracted analo- 
gously to the science spectra. All spectra were flatfielded 
using the fiber fiatfield exposures. The dedicated skyfibers 
were employed to calculate an average sky spectrum. This 
sky spectrum was then subtracted from the science spec- 
tra using the iRAF-task skytweaU* . The resulting accuracy 
(the dispersion, in a 68% confidence interval, of the medi- 
ans of the sky subtracted spectra in any exposure divided 
by the median sky) was typically better than 2%. 

Finally, the dispersion-corrected, sky-subtracted spec- 
tra were co-added to enhance the S/N using scombine, 
where the frames were weighted by their individual S/N. 
The resulting frames were then continuum-normalized or 
rectified by a polynomial fit to the regions excluding the 
Ca and other stronger absorption lines. Thus the median 
S/N ratio obtained in this way lies at 33 pixel" ^, reaching 
« 200 pixel" ^ for the brightest stars. Sample spectra of 
stars representative of our sample are shown in Fig. 3. 

2.4. Membership: Radial velocities 

In order to assess the membership of each target star, 
we made use of the individual radial velocities. The de- 
termination of accurate velocities by means of template 
spectra will be described in a forthcoming paper that will 
concentrate on the kinematic aspects of Carina (Wilkin- 
son et al., in prep.). For the moment it should suffice 
to note that the Carina dSph galaxy clearly stands out 
against the Galactic foreground contamination, at a sys- 
temic heliocentric radial velocity of ~ 220 kms"^. Mem- 
bership probabilities were derived by performing an error- 
weighted maximum likelihood fit of a Caussian velocity 
distribution between 150 km s"^ and 300 km s"^. Reject- 
ing 3(j-outliers, the fit was iterated until convergence was 
reached, resulting in a mean heliocentric radial velocity of 
223.9kms^^ and a dispersion of 7.5kms^^ (see Fig. 4). 
Taking into account our measurement uncertainties, this 
agrees very well with previous results. For example, Ma- 
teo (1998) quotes a systemic velocity of 224±3kms"^ and 
a central velocity dispersion of 6.8±1.6kms"^ (based on 

^ IRAF is distributed by the National Optical Astronomical Observator; 
In the present representation the 2a, 2.5a and 3 a limits correspond 



17 centrally located stars). Majewski et al. (2005) find a 
mean velocity of 222.8kms"^ from 61 stars that extend 
out to 1.5 nominal tidal radii. 

Depending on the chosen probability selection^" to de- 
fine a member, the total number of stars belonging to Ca- 
rina varies from 471 {2a), 501 (2.5a) to 510 (3 a). The 
number of high-probability member stars with spectra suf- 
ficiently high to be useful for abundance determinations, 
as opposed to radial velocity determinations, is somewhat 
lower, namely 409 {2a), 433 (2.5 cr) and 437 (3cr). 

Membership of red giants in the four calibration clus- 
ters was also determined via radial velocities. Only five 
stars (N4147-S370 and N4147-D57 in NGC4147; M68- 
S155 and M68-S164 in NGC4590; and M5-1046 in the 
cluster NGC 5904) were rejected due to their strongly dis- 
crepant (> 8a) velocities. 



3. CALIBRATION OF THE METALLICITY SCALE 

Determinations of acciirate spectroscopic abundances 
require a well calibrated and widely applicable reference 
scale. The infrared lines of the singly ionized calcium ion at 
A A 8498, 8542, 8662 A have become one of the spectral fea- 
tures of choice. These triplet lines are relatively easy to be 
distinguished, since they provide the strongest absorption 
lines in the near-infrared regime of RGB spectra. They 
also lie at a wavelength range where many CCDs have 
fairly high sensitivity, an advantage when trying to mea- 
sure faint sources within reasonable exposure times. How- 
ever, the CaT lines are generally considerably affected by 
strong night-sky emission lines, hampering accurate mea- 
surements for abundances in faint stars. 

The use of the CaT for measuring metallicities of red 
giants in old, metal-poor populations was first established 
for simple stellar populations such as globular clusters (Ar- 
mandroff & Zinn 1988; Armandroff & Da Costa 1991; Rut- 
ledge et al. 1997a). It was then also employed for mixed 
populations beyond the Milky Way halo, such as in dSphs 
(e.g., Suntzeff et al. 1993, Cote et al. 1999; Smecker-Hane 
et al. 1999; Guhathakurta et al. 1999; Tolstoy et al. 2001; 
Pont et al. 2004) and in irregular galaxies (e.g., Cole et al. 
2005). Based on the comprehensive catalog of CaT data 
from Rutledge et al. (1997a), Rutledge, Hesser & Stet- 
son (1997b) established a seminal calibration of Galactic 
globular cluster metallicities to the scales of Zinn & West 
(1984, hereinafter ZW) and Carretta & Gratton (1997, 
hereafter CG). Nevertheless, this empirical linear calibra- 
tion is strictly valid only for old and metal poor red gi- 
ants with [Fe/H] below — 0.5dex and ages no younger than 
10 Gyr, since this is the parameter range covered by their 
calibrating globular clusters. 

Smecker-Hane et al. (1999) carried out CaT spec- 
troscopy of 52 red giants in Carina and found that most of 
the Carina RGB stars belong to a metal poor population 
around — 2dex. Still, the occurrence of a wide range of 
ages, particularly the dominant intermediate-age popula- 
tions (Hurley-Keller et al. 1998; Rizzi et al. 2003), hence 
raises the question of an unbiased applicability of the CaT- 
technique to our Carina sample. However, recently the 
CaT calibration was extended to ages as young as 2.5 Gyr: 

to membership probabilities of p >0.0072, 0.0023 and 0.0006. 
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Based on spectroscopy of red giants in Galactic globular 
and open clusters Cole et al. (2004) have shown that the 
CaT method has very little age sensitivity and can be re- 
liably applied to composite populations with a range of 
ages. Since the dominant populations in the Carina dSph 
lie in the range of validity of the CaT method, we can trust 
that our measurements will yield consistent results for this 
galaxy. 

3.1. Equivalent widths and reduced width 

Equivalent widths of the three CaT lines were measured 
by means of the software kindly provided by A. Cole, which 
is a modified version of G. Da Costa's EWPROG-program 
(see Cole et al. 2004). We followed Suntzeff et al. (1993) in 
using line and continuum bandpasses as defined in Arman- 
droff & Zinn (1988). These were shifted to take account of 
our derived radial velocities. As was proposed by Cole et 
al. (2004), each CaT line is best fit by the sum of a Gaus- 
sian and a Lorentz profile with a common nominal line 
center, which yields a progressively better fit to the line 
wings. Theoretical models show that the wings arc most 
sensitive to changes in stellar parameters (e.g.. Smith & 
Drake 1990). The final width was then derived by sum- 
ming up the flux in the fit profile across the bandpass and 
errors were calculated automatically from the residuals of 
the fit. 

In order to perform a self-consistent analysis of our cal- 
ibration clusters and of the Carina data, we have to re- 
derive the relation between line strength and metallici- 
ties. We chose to exclude the weakest of the three lines at 
8498 A in order to minimize the errors in the final metal- 
licity measurements. Thus we define the Ca linestrength 
as 

Y:W = EW8542 + EWs662. 

The overall width T,W of the CaT lines is not only 
a strong function of metallicity, but is also affected by 

stellar cff'cctive temperature and surface gravity. Isolat- 
ing the strong sensitivity of the linestrength to metallic- 
ity is commonly achieved by defining a reduced equivalent 
width, which effectively removes the signature of surface 
gravity and effective temperature from the CaT strength 
(Rutledge et al. 1997b and references therein). Basically, 
correcting the linestrength for a star's V-magnitude over- 
comes the sensitivity to gravity, since a given star on the 
RGB can have a range of radii depending on its evolution- 
ary phase, although its mass remains nearly constant. The 
increase in radius toward the tip of the RGB results in a 
progressively lower surface gravity (and effective tempera- 
ture), where the relation between gravity and V magnitude 
is close to linear (e.g.. Cole et al. 2000). The additional 
introduction of the apparent magnitude of the horizontal 
branch, Vhb, and the subsequent measurement of the dif- 
ference between the apparent V-band magnitude of a given 
red giant and of Vh b removes any remaining dependence 
on systemic distance and reddening. 

As a comparison with the maps of Schlegel, Finkbeiner 
& Davis (1998) reveals, the Galactic foreground extinc- 
tion toward the region of Carina is rather constant (with 
an r.m.s. scatter around the mean E(B— V) of less than 
0.01 mag across the projected area of the galaxy). Thus 
differential reddening is unlikely to significantly falsify 
the use of a single HB magnitude throughout this dSph. 



In a gas- and dust-deficient dSphs like Carina the inter- 
nal reddening may be neglected. Furthermore, since the 
color range across the RGB is small and since we are us- 
ing Johnson-Cousins colors, any extinction-induced mag- 
nitude differences for stars of different surface temperature 
may be entirely ignored (Grebel & Roberts 1995). Hence 
the reduced width, defined as 

W = Y^W + p{V - Vhb), 

is basically the CaT-linestrength of an RGB star with the 
luminosity that it would have if it were located on the HB. 
This parameter is a very convenient metallicity indicator 
(Rutledge et al. 1997b), since in this form it is purely a 
function of metallicity. 

Fig. 5 shows a plot of the measured linestrengths for 
each of the calibration clusters versus the respective mag- 
nitudes above their HBs. We adopted the HB magnitudes 
of our calibrators from Ferraro et al. (1999) in order to 
place all measurements on a consistent scale. Likewise, we 
defined the HB level of Carina via the magnitude of the 
lower envelope of the observed HB distribution in the re- 
gion with 0.2< B — V <0.6. Hence we consistently use the 
determined Yhb magnitude of 20.88 mag. This is the HB 
of the old population in the Carina dSph; the core-Helium- 
burning stars of the (dominant) younger/intermediate-age 
population form the prominent red clump at ~ 20.5 (see 
Monelli et al. 2003; Smecker-Hane et al. 1994). As was 
demonstrated by Cole et al. (2004) using the HB location 
of the old HB is perfectly acceptable also for intermediate- 
age populations. We discuss the (small) expected ampli- 
tude of the resulting systematic effect in the metallicities 
of the intermediate-age stars below. 

A linear, error-weighted least-squares fit to the line mea- 
surements of all remaining 72 stars in the four globu- 
lar clusters leads to the parameterization of the reduced 
width, W', as 

W = TW + (0.55 ± 0.02) {V - Vhb), (1) 

with an r.m.s. scatter around this averaged slope of 0.18 A. 
In this equation, the given slope is the error- weighted mean 
of the individual slopes of each of the four clusters. The 
value for the slope derived by us is shallower than the 
one obtained in the original calibrations by Rutledge et 
al. (1997a), who measured it to be (0.64 ±0.02) Amag"^ 
based on the CaT line strengths in 52 Galactic globular 
clusters. The more recent analysis of Cole et al. (2004), 
who included metal-rich open clusters, obtained an even 
steeper slope of 0.73 ± 0.04. They suggest that the differ- 
ence is due to an increase of slope with the mean metal- 
licity of the calibrating clusters. Indeed, the mean metal- 
hcity of our four calibrating clusters, at —1.47 dex on the 
CG scale (-1.70 dex on the ZW scale), is ~ 0.2 - 0.4 dex 
more metal-poor than the entire sample used by Rutledge 
ct al. (1997a, b) (mean metallicities of —1.25 dex on the 
CG scale, —1.33 dex on the ZW scale, resp.), giving some 
support to the interpretation of Cole et al. (2004). Our cal- 
ibrating clusters were selected to cover the expected metal- 
licity range of the Carina dSph. Further, our analysis is 
based on only the two stronger CaT lines and uses a differ- 
ent line profile to that used by Rutledge et al. (1997a). In 
this paper, we shall consistently use our own calibration. 
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3.2. The metallicity scale 

One of our observed globular clusters, NGC 4147, was 
not included in the common scale of Carretta & Grat- 
ton (1997). In order to exploit the maximum number of 
calibrators at hand, we nevertheless retained this cluster 
in our analysis by adopting the reference mctallicities by 
Rutledge et al. (1997b), who also measured this cluster. 
These authors calibrated their sample of globular clusters 
to both the ZW and CG scales, and we also apply calibra- 
tions to both scales for the sake of consistency. Another 
calibration scale is provided by Kraft & Ivans (2003), who 
derive globular cluster metallicities based on the equiva- 
lent widths of Fe ii lines from high-resolution spectroscopy. 
Their scale yields metallicities that are about 0.2 dex more 
metal poor than those on the CG scale. However, we will 
present our results later on in terms of the calibration of 
Rutledge et al. (1997b) on the scale of CG unless stated 
otherwise, since this is the most common representation 
nowadays. 

As Fig. 6 implies our observations yield well-defined lin- 
ear relations between the reduced width W' and [Fe/H] 
on both scales. We can thus infer the metallicities of our 
target stars using the following equations, which were ob- 
tained by an error- weighted least-squares fit: 

[Fe/H]^g7_^g, = (-2.77 ± 0.06) + (0.38 ± 0.02) W^(2) 
[Fe/H]^97_^^^ = (-2.74 ± 0.06) + (0.30 ± 0.02) W[2,) 

The bottom panel of Fig. 6 shows that there is only 
little scatter around the best-fit line. The r.m.s. scatter 
was calculated to be 0.02 dex for both metallicity scales. 
The uncertainties of our metallicities were estimated by 
propagation using the measurement errors in the EWs and 
uncertainties in the c;alibration (eqs. 1 3). The mean mea- 
sured uncertainty in our metallicities amounts to 0.17 dex. 

Red giants that were observed in Carina are shown in 
the V,I]W-plane of Fig. 7, together with some isometallic- 
ity lines to guide the eye. The wide range of values in the 
diagram already indicates a large spread in metallicities. 

3.3. The possible impact of Ca abundance variations 

The "metallicity" of a galaxy is customarily represented 
by its iron abundance [Fe/H] (or [Me/H]). Strictly speak- 
ing, Ca does not necessarily trace Fe as such, and Ca and 
Fc arc formed in different nucleosynthesis processes. A 
potential inconsistency in applying the CaT technique to 
derive metallicities in a dSph in terms of [Fe/H] is the likely 
variation of the [Ca/Fe] ratio, both within the composite 
population of the dSph, and in between the calibrating 
clusters and the dSph stars. When converting from Ca to 
Fe abundances, one has to consider the different depen- 
dencies of [Fe/H] and [Ca/Fe] on the detailed star forma- 
tion histories of the galaxy. It is far from straightforward 
to determine an unambiguous metallicity scale from sin- 
gle element measurements, particularly in a galaxy with 
a strongly episodic star formation history like Carina (see 
Cole et al. 2000; Tolstoy et al. 2001 for a detailed discus- 
sion and Gilmorc & Wyse 1991 for models). 

The measured [Ca/Fe] ratios for stars in globular clus- 
ters tend to be higher than the solar value (see the com- 
pilation in Carney 1996). Indeed, the values for three of 
our calibrating clusters (each based on only a handful of 



stars; Carney 1996) are (see also Table 1): [Ca/Fe] = 
-1-0.11 ±0.04 (NGC 3201): [Ca/Fe] = +0.32 (NGC 4590; 
only two stars); [Ca/Fe] = -J-0.21±0.02 (NGC 5904). The 
few published elemental abundances for stars in the Carina 
dSph show some stars with similar values, but a tendency 
toward the lower values as expected for an extended star- 
formation history and incorporation of Fe from Type la 
supernovae. Shetrone et al. (2003) analyzed five stars that 
are known to be members from the radial velocity mea- 
surements of Mateo et al. (1993). They obtained values 
for [Ca/Fe] of: -1-0.14, -0.10, -h0.20, +0.12 and -0.02, 
with a typical uncertainty of 0.05 (see Table 4). The val- 
ues of [Fe/H] that were based on Fe lines in Shetrone et 
al.'s (2003) high-resolution spectra for these bright stars 
are: —1.60, —1.64, —1.60, —1.40 and —1.94 respectively. 

We note that uncertainties resulting from the unknown 
variations in the values of [Ca/Fe] for our target stars are 
likely to be up to 0.2 dex, based on these measured values. 
Such effects may in the future be overcome by calibrating 
W' directly onto [Ca/H] (Hosier et al. 2002), making cali- 
brations independent of star formation history, unlike the 
present technique. Our own high-resolution UVES spec- 
tra for approximately 30 stars will provide more elemental 
abundance information, which can be used to obtain an 
improved statistical estimate of the metallicity uncertainty 
(Koch et al., in prep.). 

The high-resolution measurements from Shetrone et al. 
(2003) quoted above in fact also provide an external check 
of our metallicity estimates. Comparing their results with 
our data permits us to assess the internal precision of our 
measurements. Table 4 gives a comparison of [Fe/H] val- 
ues from our analysis with the respective high-resolution 
data. The mean offset between both analyses is 0.09 dex 
(our CG scale) and 0.13 dex (ZW). After adjusting the 
individual deviations by these mean offsets, the absolute 
average difference amounts to 0.18 (0.25) dex. This offset 
is comparable to the order of uncertainty introduced by 
the [Ca/Fe] variations with respect to Galactic globular 
clusters as mentioned above, as well as our random mea- 
surement errors. Since we here juxtapose our Ca-based 
metallicities with the true iron abundances from Shetrone 
et al. (2003), the above quoted average difference well re- 
flects the calibration uncertainty. 

3.4. Additional potential sources of uncertainty 

Furthermore, while ideally one would like to associate 
each RGB star with the appropriate HB level correspond- 
ing to a population of the age and metallicity of that RGB 
star, in practice we do not have sufficient information. Our 
adoption of a single HB for all the stars has the advan- 
tage of providing us with a well-defined magnitude and 
allows us to estimate the associated uncertainties. Since 
we adopted the HB appropriate for the oldest stars, and 
hence the least luminous HB, there will be a systematic 
effect. The expected result from normalizing with an HB 
of apparent magnitude VHB,true + ^Vrb, with VHB,true 
the appropriate HB level for a given star, can be seen by 
the appropriate substitution in the relations for reduced 
width and [Fe/H] (eqs. 1,3). We find a resulting error in 
[Fe/H] of -0.55 x 0.38 x AVhb dex. As noted earlier, the 
core He-burning phase of the intermediate-age population 
in the Carina dSph forms a prominent red clump, some 
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0.4 mag brighter than the old HB used in the normaliza- 
tion; thus for these stars AVhb ^ 0.4. The resulting bias 
in [Fe/H] is —0.07 dex, i.e., we estimate that these stars 
may be —0.07 dex too metal-poor. 

The small amplitude of this effect is confirmed by the 
analysis of Cole et al. (2004), where they experiment by 
assigning HB levels (absolute magnitudes) randomly to 
stars in different clusters. The variation in the luminos- 
ity of the clusters' HB levels is 0.4 mag (sec their Ta- 
bic 1). Using their derived slopes for the calibration re- 
lations, we find that the total amplitude of the change in 
metallicity estimates from using the appropriate HB level 
is 0.1 dex. This is just what is seen in their Figure 6. 
Basically the variation of HB level with age and metallic- 
ity amounts to ~ 0.1 dex in metallicity uncertainty, and 
in our case for most of the stars results in an underesti- 
mate by this amount. Of course we cannot say for which 
stars we arc underestimating the metallicity without hav- 
ing age information. As we discuss below, in principle, 
with robust accurate models of the upper RGB, one could 
iteratively solve for age and metallicity. We did not at- 
tempt this procedure since the maximum amplitude of the 
effect, 0.1 dex, is within the uncertainties of the models. 

Of the target stars in the central field, 46 were observed 
repeatedly, i.e., during all of the runs. In these cases the 
finally quoted metallicity is the value derived from a com- 
bined spectrum, including all of the respective single spec- 
tra. However, in order to get an estimate of the random er- 
rors of the measurements, we compare in Fig. 8 the metal- 
licities determined from spectra obtained on the individual 
observing runs in 2003 February and March to metallici- 
ties derived from combined spectra from all the available 
observations (see also Table 5). Some of the entries have 
metallicity estimates from the two runs differing by sev- 
eral sigma. It turns out that the median difference ([Fc/H] 
(1st run) — [Fe/H] (2nd run)) between individual metal- 
licity estimates amounts to 0.08 dex and the r.m.s. scatter 
of the different estimates from both runs totals 0.57 dex. 
Rejecting those objects deviating more than two standard 
deviations, namely two targets with particularly low S/N, 
from the sample we find a respective r.m.s. of 0.3 dex. 
This is larger than our overall random error estimate of 
0.17 dex. But as noted before, this is mainly due to the 
first run being plagued by poor weather conditions, which 
was taken into account by weighting spectra by their S/N 
ratio during the c;o-addition. Consequently, the mean un- 
certainty of 0.24 dex for the targets observed during the 
first run compares to 0.13 dex for the same objects tar- 
geted in the second run. 

4. THE METALLICITY DISTRIBUTION 

Metallicities derived for each individual star are given in 
Table 6. The resultant MDFs (see Fig. 9) turn out to be 
almost indistinguishable when using different membership 
selections (i.e., cutting at different a levels in the radial 
velocity distributions), except for the varying total num- 
ber. We used both a K-S test and the Kuiper statistics 
(Press et al. 1992) to test how different the MDFs are for 
different selection criteria. The Kuiper statistics is more 
sensitive to differences in the MDFs' tails. The MDFs re- 
sulting from a la, 2a, and 3cr-cut are consistent with each 
other at a level of greater than 99.9%. Hence we chose to 



maintain the limit of 3 ct as our membership criterion. 

4.1. The width of the MDF 

Histograms of the resulting metallicities on each of the 
two calibration scales for the sample of the 437 Carina red 
giants are displayed in Fig. 9. These show a broad distri- 
bution, peaked at a mean metallicity of —1.72 ± 0.01 dex 
(CG scale) or -1.90 ± 0.01 (ZW scale), respectively. The 
peak of the MDF is, however, more pronounced for the 
case of the CG-scale. The entire distribution's formal full 
width at half maximum is 0.92 dex. The full metallicity 
range is ^ 3.0 dex, as can be seen in Fig. 9. This spread 
can in part be attributed to the measurement uncertain- 
ties, but also an occurrence of several subpopulations with 
different peak metallicities cannot be rejected as a further 
source of the broadening (sec below) . 

The formal 1 cr-width of the distribution is 0.39 dex, but 
since a best-fit single Gaussian is simply a mathematically 
convenient description, some of this width will also be due 
to measurement and calibration uncertainties. Subtract- 
ing our best estimates of these, namely 0.17dex for the 
measurement error and 0.02 uncertainty in the slope of 
the calibration in quadrature, we estimate the width of 
our distribution to be a =0.33 dex. One should keep in 
mind, however, that the combined effects of calibration un- 
certainties from the [Fe/H] calibration via reference clus- 
ters through systematic effects due to varying [Ca/Fe] and 
normalizing the HB level may amount up to 0.22 dex in to- 
tal. Accounting also for these systematic errors, we obtain 
an estimate of the width of the MDF of ct = 0.25 dex. 
Furthermore, we recall that the dominant population in 
Carina is of intermediate age, reducing the possibility of 
star-to-star random variations of uncertainties introduced 
by the chosen HB level. Thus the width of the distribution 
is in fact dominated by physical variations in metallicity. 

In the extreme tails of the MDF we find stars with 
metallicities approaching —3 dex and near-solar metallic- 
ity, respectively. Note that our calibrating globular clus- 
ters only cover the initially expected range of metallicities 
in Carina, i.e., metallicities from ~ —2 to ~ — 1 dex (Ta- 
ble 1). Hence the metallicities for the stars in the tails 
of the MDF are based in extrapolation of our calibration. 
Follow-up spectroscopy of these stars would be desirable 
in order to check whether this extrapolation yields reliable 
results and to further elucidate the chemical properties of 
these stars. If the overall [Fe/H] values for these stars are 
correct, they would be as metal-poor as the most metal- 
poor red giants found in other nearby dwarf galaxies (e.g., 
Shetrone et al. 2001). At the metal-rich end, a handful of 
them would be as metal-rich as the metal-rich population 
in the Sgr dSph (Bonifacio et al. 2004). 

4.2. Comparison with previous studies 

We can now compare the results of our large sample with 
previous results for large samples of red giants in Carina. 
From a sample of 52 RGB stars using the same method 
as here, Smecker-Hane et al. (1999) found the mean spec- 
troscopic metallicity of Carina to be —1.99 ± 0.08 dex. 
These authors have not published individual metallicities, 
so we cannot perform a direct comparison between stars 
in common to both samples. Rizzi et al. (2003) analyzed 
the color distribution across the RGB, including a sta- 
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tistical correction for the star-formation history inferred 
from the main sequence turn-off region. They derived a 
photometry-based average metalhcity of —1.91 ± 0.18 dex. 
These results on Carina are in reasonable agreement with 
our data if the quoted uncertainties and the widths of the 
distributions are taken into account. 

Yet, metallicity estimates via purely photometrical 
isochrone fitting of RGB stars without age information 
lead to mean values that are typically more metal-rich by 
approximately 0.2 dex, as, e.g., in MoncUi et al. (2003). 
These authors obtain their result by a direct comparison 
with the Galactic globular cluster NGC 1904, a single- 
age, single-metallicity system. For this object, different 
values for [Fc/H] are found in the literature, reaching 
from —1.37 dex (Rutledge et al. 1997b on the CG scale) 
to — 1.69dex (loc. cit., ZW-scale), whereas Harris (1996) 
cites an average metallicity of — 1.57dex from all available 
sources. Still, one has to keep in mind that such a purely 
photometric fit to the RGB, without any input from an 
age distribution and age-metallicity relation, is subject to 
large uncertainties, as is underscored by Fig. 10. An addi- 
tional source of uncertainty is the ability (or lack thereof) 
of many isochrone models to correctly reproduce the slopes 
of the red giant branch for very low or high metallicities 
(e.g., Fig. 5 in Grebel 1997 and Fig. 5 in Grebel 1999). 

4.3. The shape of the MDF 

The complex morphology of the CMD of the Carina 
dSph - both the multiple turn-off regions and the different 
He-burning phases - is compelling evidence for strongly 
varying star formation rates, in fact even for distinct 
episodes. The metallicity distribution of low-mass stars, 
as derived here, is rather insensitive to the details of the 
star formation history, particularly if elements dominated 
by massive stars are used (cf. Pagel 1997, Section 8.3.7). 
However, we have calibrated our CaT metallicities onto 
iron (assuming the same [Ca/Fe] in the calibrators as in 
our stars), which in systems with extended star formation 
can have a significant contribution from long-lived stars 
through Type la supernovae. Thus one can have enrich- 
ment in Fe without accompanying star formation. Hence 
with a star-formation history that consists of bursts of star 
formation with a real hiatus between the bursts, the Fe 
distribution can have "gaps" (no stars formed while Type 
la supernovae enriched the gas, prior to the next burst; 
cf. Gilmore & Wyse 1991). In order to test whether any 
substructure reflected in the MDF may provide evidence 
of star formation events, we applied the Shapiro- Wilk test 
for normality of a distribution (e.g., Royston 1995) to our 
observed binned MDF. As a result we can reject the hy- 
pothesis of an underlying single Gaussian normal distribu- 
tion with a significance of 96%. On the other hand, the 
observed MDF is not expected to be strictly Gaussian, 
since star formation events do not naturally produce such 
normal distributions. In addition, we applied a KMM test 
for multiple populations imprinted in the MDF (Ashman, 
Bird, & Zepf 1994) to our unbinned data. Again, given our 
measurement uncertainties and the non-normality of as- 
trophysical star formation events, one has to keep in mind 
that such a decomposition into multiple single Gaussian 
sub-populations is purely a formal procedure at this stage 
that does not necessarily have a physical meaning. We find 



that the case of four populations is significantly preferred 
against a one-population model at the 98.1% confidence 
level. Such a potential multi-population star formation 
history could also easily account for the MDF's metal-rich 
extension. Yet a proper interpretation of the observed 
metallicities requires detailled modeling of the chemical 
evolution to which we will return in Sect. 6 below. 

4.4. An age-metallicity relation? 

The (lack of a) relation between the color of our tar- 
gets and their metallicity is demonstrated in Fig. 10. This 
figure shows that there is no obvious trend in the sense 
that predominantly metal-poor stars are not concentrated 
along the blue ridgeline of the RGB, and the metal-rich 
stars are not concentrated at the red ridgeline (in contrast 
to the LMC data of Cole et al. 2005, where such a trend is 
clearly visible) . A metallicity-color trend would have been 
expected to result from a direct comparison with fiducials 
of globular clusters of known metallicities, i.e., for popu- 
lations of roughly the same age but diff'erent metallicities, 
as indicated by the RGB fiducials for four globular clus- 
ters of Sarajedini & Layden (1997) in Fig. 10 (left panel). 
The ridgelines from Sarajedini & Layden cover a wide 
range of colors in the CMD, corresponding to a metal- 
licity range from —2.02 to —0.78 dex for old and nearly 
coeval populations. On the right-hand side we show three 
isochrones with ages and metallicities that roughly corre- 
spond to what may be expected from the three main star 
formation episodes in Carina. Again there is no clear cor- 
relation with the data points. Hence it is clear that one 
cannot derive the metallicity of an individual star from 
its color and magnitude on the RGB when dealing with 
a mixed-age population as present in Carina, nor can a 
metallicity distribution directly be inferred from the color 
distribution in this case. 

The existence of an age-metallicity relation can be in- 
ferred from the fact that the giants appear to be equally 
widely spread across the RGB regardless of their metallic- 
ity. At a given age, metal-poor stars are in general bluer 
than their metal-rich counterparts owing to the reduced 
opacity in their photospheres. Conversely, the older stars 
tend to have redder colors at a given [Fc/H]. Hence, in or- 
der to produce the counteracting trend with both metal- 
rich and metal-poor stars at the same locus in the CMD, 
the metal-poor stars necessarily have to be older, whereas 
the more metal-rich giants have to be of younger age. 

If one had reliable models of the upper RGB in the 
observational plane, one could in principle derive age es- 
timates for our stars of known metallicity, from finding 
the best-matched isochrone at that metallicity. However, 
the small offset in color with age at the relevant inter- 
mediate ages, combined with sensitivity to the largely un- 
known elemental abundance mix and He content limits the 
usefulness of this exercise. An age distribution and age- 
metallicity relations are better derived from turn-off stars, 
provided that the photometry is at a sufficient level of ac- 
curacy. Such an analysis with high-resolution spectra will 
be subject of a forthcoming paper (Koch et al., in prep.). 

5. RADIAL VARIATIONS 

Previous analyses based on wide-area photometric sur- 
veys and thus derived CMDs have shown that the stellar 
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populations in the Carina dSph feature different spatial 
distributions. Mighell (1997) proposed a strong central 
concentration of intermediate-age stars from a compari- 
son between a deep HST CMD of the center and ground- 
based data that covered a wider area. This was confirmed 
by Harbeck et al. (2001) via wide-field photometry. Har- 
beck et al. found a pronounced central concentration of 
the prominent intermediate-age red clump stars in com- 
parison to the old HB. This was subsequently also seen in 
the photometric study of Monclli et al. (2003). Harbeck 
et al. (2001) discarded the presence of a significant radial 
gradient in Carina's old HB morphology out to its nominal 
tidal radius. 

Our large dataset allows us to address the issue of pop- 
ulation gradients by looking for radial gradients in the 
MDFs. For this purpose we investigated the MDFs as 
a function of elliptical radius r = [x'^ + y^/(l — e)^]^/^, 
where x and y denote the distances along the major and 
minor axis, respectively^^. There is no obvious strong ra- 
dial gradient in the resulting scatter plot of metallicity 
against elliptical radial coordinate, as shown in Fig. 11 
(top panel). A simple linear fit yields a gradient of 
(-0.0047 ± 0.0008) dexarcmin"!. 

However, as Fig. 11 (bottom left) implies, there is a 
clear radial tendency in the respective MDFs to be seen 
when the data are split into three radial bins. The an- 
nuli have been chosen to cover everything within one core 
radius {ic = 8.'8), then everything within Tc to 2rc, and 
from 2rc to the nominal tidal radius at 28.'8. The ratio of 
the total number of stars in the annuli is approximately 
5:4:1. Fig. 11 shows the single MDFs, represented by nor- 
malized generalized histograms, which were obtained from 
the density distributions of our [Fe/H] measurements by 
convolution with a Gaussian distribution accounting for 
the individual measurement errors. 

The MDF of the innermost region is already rather 
broad with a suggestion of both lower and higher metallic- 
ities and an indication of a secondary population around 
approximately —2 dex. However, the MDF of the middle 
annulus is peaked at a metallicity that is at least 0.1 dex 
lower than in the inner region and shows an increase of 
the potential secondary peak at ~ —2 dex. Considering 
the overall uncertainties in the metallicity measurements, 
a shift by, e.g., 0.1 dex may not appear significant. How- 
ever, if we take into account that the uncertainties of the 
metallicity measurements should be the same everywhere 
in Carina and that there is no obvious reason why these 
uncertainties should introduce systematic offsets from one 
bin to another, these differences may be meaningful. 

A K-S test revealed a probabihty of 98.2% that the MDF 
from the innermost ring is drawn from the same parent 
population as the entire sample. Based on the same sta- 
tistical test, we cannot reject the hypothesis that also the 
middle annulus is drawn from the same population as the 
innermost annulus and the entire sample, resp. (at signif- 
icances of 35% and 50%). 

The MDF of the outermost region, on the other hand, 
shows several peaks, the broadest and most distinct of 
which is more metal-poor by ^0.3 dex than that the one 
seen in the innermost region. In contrast, the location of 
a second peak at an [Fe/H] of —1.6 dex is consistent with 

Here we adopted a position angle of (65 ± 5)° and an ellipticity of e 



the MDFs' peaks in the inner parts of the galaxy to within 
the measurement uncertainties. Using a K-S test we find 
that the probability that the outer and inner MDFs are 
drawn from the same distribution is less than 4% and we 
can therefore reject this hypothesis. Moreover, the proba- 
bility of the outer region being compatible with the entire 
sample is at the 6% level. Hence there do seem to be indi- 
cations of spatial (and indeed radial) variations in Carina's 
stellar populations. 

This gradient is also reflected in the cumulative metal- 
licity distributions shown in the bottom right panel of 
Fig. 11. In a complementary manner, a cumulative plot of 
the giants' galactocentric radii clearly indicates a prefer- 
ential central concentration of the metal-rich component 
(Fig. 12) compared to the more metal-poor population - a 
K-S yields a probability that the metal-poor and metal- 
rich populations take the same spatial distributions of 
10%. 

The evidence of a concentration of the metal-rich com- 
ponent in Carina that we found is consistent with ear- 
lier findings that the intermediate-age stars favor the in- 
ner regions, in particular if one takes into consideration 
that we have likely underestimated their metallicity by 
~ 0.1 dex due to the adoption of the old HB in the re- 
duced CaT width. Thus, if the metal-rich component can 
be attributed to the intermediate-age population and, vice 
versa, the metal-poor RGB stars in Carina are to be identi- 
fied as predominantly belonging to an old population, then 
the radial variations and progressive shifts of the peaks 
found in the spectroscopically determined MDFs under- 
score previous findings of centrally concentrated starbursts 
based on CMD analyses (Harbeck et al. 2001). A reason 
for such a concentration could be stronger cooling and dis- 
sipation of more metal-enriched gas. It may also indicate 
that star-forming material can be more easily retained at 
the center of a dSph's shallow potential well. 

6. IMPLICATIONS FOR CHEMICAL EVOLUTION 

As the present-day observed MDF of a stellar system 
reflects its entire integrated history it is relatively insensi- 
tive to the star-formation history (SFH). The SFH cannot 

be unequivocally derived unless these data are combined 
with accurate ages (and preferably also individual elemen- 
tal abundances to constrain the actual modes of star for- 
mation, see, e.g., Shetrone, Cote, & Sargent 2001). How- 
ever, the gross properties of the MDF can provide insight 
into certain aspects of the chemical evolution of a sys- 
tem, such as the importance of gas flows. Following long- 
standing practice, we gain a basic insight by consideration 
of variants of the simple model of chemical evolution (e.g., 
Pagel 1997, Chapter 8) and refer the detailed stochastical 
modeling to a future paper (Wyse et al., in preparation). 

In such basic models, the mean metallicity of long-lived 
stars at late stages (close to gas exhaustion) is represen- 
tative of the true nucleosynthctic yield, p, if all metals 
are retained in the system. If the system is not closed 
and if outflows are possible, the mean metallicity reflects 
the 'effective yield", a combination of true yield and ef- 
fects of loss of metals (see Hartwick 1976): outflows at a 
constant factor c times the star-formation rate reduce the 
mean metallicity below the true yield by approximately 

= 0.33 ± 0.05 (Irwin & Hatzidimitriou 1995). 



10 



Koch et al. 



this factor c, PeS ~ p/c. Thus, assuming an invariant stel- 
lar initial mass function (and thus invariant true yield p) , 
the low peak metallicity of our derived MDF, compared to 
say the G dwarfs in the solar neighborhood, which peaks 
at ^ —0.2 dex (e.g., Nordstrom et al. 2004 and references 
therein) suggests a strong influence of outflows, presum- 
ably driven by Type II supernovae. Indeed, this has been 
the basic model for the evolution of gas-poor dwarf galax- 
ies since Sandage (1965). 

6.1. Outflows 

A modified simple model with flows overpredicts the 
number of metal-poor stars (leading to a "G-dwarf prob- 
lem" just as in applications of this model to the solar 
neighborhood). Having normalized the fit of our mod- 
els in Fig. 13 to the total number of stars observed, 
the discrepancies between the prediction and observations 
are marginal at —2.5 dex, but become significant below 
—2.7 dex, where the deviations progressively exceed the 
2 cr-level. The fit may be improved by adopting some non- 
zero initial metallicity (or prompt initial enrichment (PIE), 
e.g., Tinsley 1975), understood in the present-day context 
as being from pre-galactic Population III stars. The pre- 
dicted metallicity distributions for a simple model with 
outflows, and with PIE, are shown compared to the full 
Carina data set in Fig. 13. The respective best- fit pa- 
rameters to our data are peff = 0.028^0 for the simple 
model, and Pe// = 0.029 Zq and Zinu = 0.0013 in 
the case of PIE. One should note that these models retain 
the instantaneous recycling approximation and so are not 
expected to provide a perfect fit to our iron-based MDF. 
Further, the radial variations discussed in Section 5 im- 
ply that the system is not in fact well- mixed at all times, 
another assumption inherent to the simple model. Given 
these caveats, the models provide an adequate fit. There 
are still some discrepancies present, such as the prominent 
main peak at [Fc/HJcg ^ —1.7 dex, which deviates by 
more than 2a (taken as ^/N) from the model curves and a 
single Gaussian. Note, however that this peak is not dis- 
tinct in the MDF on the ZW scale. As noted above, any 
potential substructure may reflect variations in star for- 
mation activity. That such variations must have occurred 
is obvious from the presence of multiple main-sequence 
turn-offs and the decline in star formation activity in be- 
tween. Indeed, not just outflows are implicated, but also 
inflows to sustain star formation over many Gyr (see, e.g.. 
Silk, Wysc & Shields 1987, Lin & Murray 1998; Carigi, 
Hernandez & Gilmore 2002; Dong et al. 2003). 

6.2. Infall and the gravitational potential 

Mac Low & Ferrara (1999) investigated starburst-driven 
mass loss from dwarf galaxies. They distinguish between 
gas that is blown out but retained by the gravitational po- 
tential well of their host galaxy, and gas that is entirely 
blown away and lost from the galaxy. Lin & Murray (1998) 
and Dong et al. (2003) present a model for episodic star 
formation in dwarf galaxies like Carina according to which 
enriched gas is blown out but eventually re-accreted once 
it has cooled down sufficiently, leading to a new episode 
of star formation. Smecker-Hane et al. (1996) argued that 
episodic accretion of fresh gas seemed unlikely since the 
mass of Carina's dark halo seemed to be too low; also Ca- 



rigi et al. (2002) note that in the case of external accretion 
the previously measured tidal radius (Irwin & Hatzidim- 
itriou 1995) needs to be larger, i.e., the potential well needs 
to be deeper to permit the accretion of external gas. 

Odenkirchen et al. (2001) demonstrated that tidal radii 
inferred from photographic data may underestimate the 
angular extent and nominal limiting radius of dwarf galax- 
ies. For Carina we may expect a firmer limit on the dark 
matter halo from the analysis of its radial velocity dis- 
persion profile (Wilkinson et al., in prep.). Here we note 
that our spectroscopic data indicate that four red giants 
whose radial velocities are consistent with membership in 
Carina lie beyond of its nominal limiting stellar radius 
(Fig. 11, top panel). Three of these stars have metallici- 
ties consistent with the trend seen in the outskirts of Ca- 
rina. Claims suggesting the existence of "extratidal" stars 
around Carina have also been made by Kuhn, Smith, & 
Hawley (1996) and Majcwski et al. (2000, 2005). Whether 
this implies that Carina is not in virial equilibrium and 
in the process of losing/having lost substantial amounts of 
mass as suggested in these studies is beyond the scope of 
the present paper. However, if Carina were a remnant of 
a once substantially more massive galaxy then the simple 
chemical models discussed here would need to be modified. 

6.3. A model with outflows and infall 

A more sophisticated model (Lanfranchi & Matteucci 
2004) is also plotted in Fig. 13. This special model uses as 
input a SFH consistent with the CMD and hence incorpo- 
rates two major epochs of outflows and respective winds 
associated with two main star formation episodes, each 
with a duration of 3 Gyr. The crucial model parameters 
were fit such as to reproduce Carina's total stellar mass, 
the gas content and observed [a/Fe] abundance ratio for 
five stars. The distribution resulting from this model holds 
strictly for the true elemental iron abundance as mea- 
sured from high-resolution spectroscopy, rather than over- 
all metallicity that we derived from our CaT measurements 
(but see the discussion in Section 3.3). The Lanfranchi & 
Matteucci model for Carina shows one dominant metallic- 
ity peak, demonstrating the insensitivity of the resulting 
MDF to the details of the SFH, even for Fe. There is 
again acceptable agreement between this particular model 
(in a non-instantaneous recycling approximation) and our 
observations, particularly when observational uncertain- 
ties arc included (right panel, Fig. 13). However, as be- 
fore there are also remaining discrepancies, e.g., associated 
with the metal-rich and metal-poor tails (with a persisting 
G-dwarf problem). 

The Carina model of Lanfranchi & Matteucci tries to re- 
produce individual element abundance measurements for 
five stars and has hence constraints only over the narrow 
metallicity range of these five stars. Lanfranchi & Mat- 
teucci's "best" model (see Fig. 13) requires a wind effi- 
ciency of some seven times the SFR and high star forma- 
tion efficiency, leading to solar metallicities within only 
a few Gyr. This is at odds with our MDF in combina- 
tion with the photometrically inferred ages of the different 
star formation episodes in Carina (e.g., Hurley-Keller et 
al. 1998; Monelh et al. 2003). Lanfranchi & Matteucci 
point out that their results also differ from the SFH de- 
rived by Tolstoy et al. (2003) on the basis of the five stars 



Metallicities in Carina 



11 



mentioned before. Adding more data points (Koch et al., 
in prep.) should aid in imposing stronger constraints. 

The spatial variation of the SFH of the Carina dSph 
and other dSphs (e.g., Harbeck et al. 2001), plus evidence 
for cold sub-structures in dSphs (e.g., Klcyna et al. 2004; 
Wilkinson et al. 2004) suggest that star formation and 
chemical evolution should be modeled on small scales. We 
will present the details of such models in a future paper 
(Wyse et al., in prep.). 

7. SUMMARY AND DISCUSSION 

Carina is the only dSph known to have undergone dis- 
tinct, well-separated episodic star formation. The domi- 
nant episodes of star formation took place at intermediate 
ages. It is not yet understood what caused the repeated 
cessation and delayed re-start of star formation in this 
enigmatic dSph. In the framework of a VLT Large Pro- 
gramme aiming exploring the chemical evolution, kinemat- 
ics, and dark matter content of Galactic dSphs, we have 
derived CaT metallicities for 437 red giants in the Galac- 
tic dSph Carina. Our spectroscopic sample exceeds the 
previously largest published dataset (Smecker-Hane et al. 
1999) by a factor of ~ 8. Our targets cover the entire sur- 
face area of Carina and allow us to investigate also radial 
metallicity variations. 

While we are using Ca as a tracer of overall metallicity 
("[Fc/H]"), we are wcll-awarc of the different dependen- 
cies of [Fe/H] and [Ca/Fe] on the detailed star-formation 
history of a galaxy, and of possible differences between Ca 
and Fe abundances in globular clusters and dwarf galaxies. 
We compare our results to direct Fe measurements from 
the few existing high-resolution spectra in the literature 
and arrive at the same conclusions as previous authors: 
The widely used CaT method is indeed a powerful tool 
for efficient measurements of overall metallicities and for 
studies of overall trends in the metallicity distribution of 
intermediate-age and old populations in nearby galaxies as 
presented here. However, uncertainties in the metallicity 
measurements for individual red giants may amount to up 
to more than 0.2 dex when using the CaT method. 

The resulting metallicity distribution function shows a 
wide range of metallicities spanning ~ 3.0 dex in [Fc/H]. 
The MDF peaks at a mean metallicity of —1.72 ± 0.01 dex 
(CG scale) or -1.90 ± 0.01 (ZW scale), respectively. Our 
mean metallicity on the CG scale is slightly more metal- 
rich as compared to earlier spectroscopic CaT studies (e.g., 
Smecker-Hane ct al. 1999). The entire distribution's for- 
mal full width at half maximum is 0.92 dex. We proposed 
that some of this spread might be attributed to the oc- 
currence of an extended SF and self-enrichment, resulting 
in several subpopulations with different peak metallicities. 
The question whether these subpopulations of different 
metallicities are manifestations of corresponding distinct 
episodes of SF requires a combination of deep photomet- 
ric and spectroscopic data and detailled modelling and will 
be left for a future paper. The MDF and its peaks change 
as one goes to larger distances from Carina's center. Gen- 
erally, it shows a trend for more pronounced metal-poor 
subpopulations with increasing distance. 

We found a higher concentration of metal-rich, presum- 
ably intermediate-age stars toward the inner regions of 
the galaxy as compared to more metal-poor stars. Metal- 



poor stars are, however, detected throughout the entire 
galaxy. This trend matches the observed central concen- 
tration of intermediate-age stars in Carina (Harbeck et al. 
2001), supporting that the more metal-rich stars are also 
the younger ones. Another argument in support of an age- 
metallicity relation in Carina is the location of the more 
metal-rich stars on the RGB: Their metallicities in combi- 
nation with their photometric colors are inconsistent with 
these stars belonging to the old population. 

Our data thus confirm earlier suggestions that Carina 
shows an age-metallicity degeneracy in the sense that 
higher metallicities counteract the effects of younger ages. 
This conspiracy in turn leads to the observed narrow RGB. 
These findings and the considerable metallicity spread rule 
out scenarios in which Carina's unusual episodic star for- 
mation history was caused by the repeated accretion of 
pristine, unenriched gas. 

Comparing our MDFs to simple chemical evolution 
models such as a closed-box model and prompt initial en- 
richment reveals a G-dwarf problem at low metallicities, 
but a fairly close approximation of the shape of the MDF 
at high metallicities. The metallicity peak of the MDF it- 
self is not reproduced by these simple models. One should 
bear in mind that the peak is most prominent on the CG 
scale, but does not show up in the MDF based on the scale 
of ZW. The early star formation in Carina that can still be 
traced by present-day red giants must have occurred from 
pre-cnriched gas (e.g., by Population III stars). A better 
match to the MDF can be obtained through more sophis- 
ticated models that include infall and outflows. The steep 
risc^ toward higher metalliciti(^s and its culmination in a 
prominent peak are qualitatively well reproduced by the 
models of Lanfranchi & Matteucci (2004), but the subse- 
quent decline toward even higher metallicities is too rapid 
compared to our observationally derived MDF. Still, the 
qualitatively good overall agreement between these mod- 
els and the observations is encouraging. One possible ex- 
planation for the repeated onset and cessation of the ex- 
tended star formation episodes in Carina is the possible 
re-accretion of previously blown-out, but not blown-away, 
enriched gas (sec, e.g., Lin & Murray 1998; Dong et al. 
2003). The predicted time scales for blow-out, cooling, and 
re-accretion are consistent with the observed long pauses 
in star formation activity. 

Refined chemical evolution models will be presented 
in a forthcoming paper (Wyse ct al., in prep.), comple- 
mented by results from high-resolution spectroscopy and 
a new derivation of the star formation history of Carina 
that combines spectroscopic and photometric information 
(Koch et al., in prep.). This will also lead to a quantita- 
tive age-metallicity relation. Furthermore, future papers 
will investigate whether different subpopulations in Carina 
show different kinematics (Wilkinson et al., in prep.). 

A handful of the red giants whose radial velocities 
make them likely Carina members lie beyond the nominal 
tidal radius of Carina. Spectroscopic evidence for poten- 
tial "extratidal" stars has also been found by Majewski 
et al. (2000, 2005), raising the question whether Carina 
is a dark-matter-dominated, boimd dSph with an even 
larger limiting radius than previously thought (see, e.g., 
Odenkirchen et al. 2001) or a disrupted remnant. Carina's 
orbit as inferred from its proper motion indicates that 
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its perigalactic distance to the Milky Way is rather close 
(~ 20 kpc), suggesting that Carina needs to be strongly 
dark-matter dominated in order to have survived these 
passages (Piatek et al. 2003). The radial velocity disper- 
sion profiles that will be derived as part of of our VLT 
Large Programme will help to elucidate the extent of Ca- 
rina's dark matter halo (Wilkinson et al., in prep.). 

In spite of its unusual star formation history, Carina 
fits in well with the mctallicity-luminosity relation of 
other dSph galaxies. The mctallicity-luminosity relation 
of dSphs is offset from that of dirrs such that dSphs have 
higher mean stellar metallicities for a fixed optical lumi- 
nosity. This offset suggests that star formation in young 
dSphs was more vigorous than in young dIrrs (Grebel et 
al. 2003). As is the case for other dSphs, Carina is thus 
unlikely to simply be a stripped former irregular galaxy. 
However, Piatek et al. suggest that Carina is currently 
near apogalacticon and that its pericentric passages bring 
it very close to the Milky Way, while its star formation his- 
tory docs not seem to be correlated with its orbital motion. 
Nonetheless, Carina's orbit would seem to suggest that in 
addition to intrinsic effects such as feedback, its evolu- 
tion must also have been influenced by externally induced 
mass loss. Such mass loss may have been induced through 



tidal interactions with the Milky Way and/or through ram 
pressure stripping when passing through the denser re- 
gions of the Galaxy's gaseous halo. There appears to be 
a correlation between present-day Galactocentric distance 
and the fraction of intermediate-age populations (van den 
Bergh 1994; Grebel 1997, 1999). Carina with its distance 
of 94 kpc and its prominent intermediate-age population 
seems to fit into this correlation rather well. However, one 
may ask whether this is merely fortuitous for this galaxy,if 
one considers Carina's proposed short pericentric distance 
and the resulting eccentric orbit. Clearly, a comprehen- 
sive exploration of Carina's history will have to take into 
account both internal and external influences, its detailed 
chemical evolution, its age structure, its mass and mass 
distribution, its internal kinematics and its orbit. 
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port by the Swiss National Science Foundation through 
grant 200021-101924 and 200020-105260. MIW acknowl- 
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Flo. 1. — Region of the red giant branch in Carina, from which our targets were selected. The photometry shown here is EIS photometry 
transformed into the Landolt-Stetson UBVRI system. The soUd lines delineate the color cuts imposed to select the targets and open squares 
denote the confirmed radial velocity members of Carina. 
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Fig. 2. — Location of our targets (shown as open squares) relative to Carina's geometric centre {ao,5o), superposed on a stellar map of the 
galaxy and its surroundings. The circles circumscribe the field of view of the FLAMES spectrograph of 25'. The ellipse indicates the nominal 
tidal radius, ellipticity, and position angle of Carina as inferred by Irwin & Hatzidimitriou (1995). 
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Fig. 3. — Sample spectra of Carina red giants centered on the near-infrared Call triplet region. The three spectra are representative of the 
observed range in metallicities. 
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Fig. 4.— Left 
of ~223kms-l. 



panel: Distribution of radial velocities for all stars in our sample. Right panel: The region around Carina's systemic velocity 
Also indicated is the best fit of a single Gaussian. 
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HB 



Fig. 5. — Equivalent widths of red giant branch stars in our four calibrating globular clusters jjlottcd against their magnitude above their 
respective horizontal branches. Shown as a solid line is the best-fitting line for each cluster, based on a common slope of 0.55Amag~^. 
The symbols denote NGC 5904 (open circles), NGC320f (open squares), NGC4f47 (open triangles) and NGC 1.590 (open diamonds). A 
representative Icr errorbar is indicated at the lower left. The metallicities of the globular clusters decrease from top to bottom. See Table 1 
for their [Fe/H] values on different metallicity scales. 
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Fig. 6. — The top panel displays the inctallicity of the four calibration clusters against their reduced width, both calibrated against the 
reference values of Rutledge et al. (1997b) on the scale of Zinn & West (stars) and Carretta & Gratton (triangles). The residuals of the linear 
best fit (dashed lines) are plotted in the bottom panel. 
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Fig. 7. — Distribution of the Carina targets in the V,SW-plane. As the lines of constant metallicity indicate, the galaxy exhibits a larg 
spread in [Fe/H]. A typical errorbax is shown left. 



20 



Koch et al. 



-0.8- 



-1 



-1.2 



J -1.4- 



o 



o -1.6h 
o 



-1.8- 



-2 - 



-2.2 - 



-2.4 



(D 
O 



1 - 







I — e — I 



I — o 



I e 1 



I — o 



/ 



/ 



/ 



/ 



I — — I 



Q) — I 



^5 



/ -- 



-I— I 



W0 



/ 



/ 



/ 



/ 



r 



/ 



/ 

/\ 1 



-M O 1 



^) 1 



1 



(> 




() 



() 



-2.5 



-2 -1.5 

[Fe/H] (Combined) 



-1 



Fig. 8. — Comparison of the stars observed during different observing runs: Shown arc metallicitics derived from spectra of the first run 
(labeled carlj) versus [Fe/H] as measured from the combined spectra of all runs. The dashed line is unity and the bottom panel displays the 
respective residuals. 
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100 - 




Fig. 9. — Histograms for metallicities of the RGB stars in Carina. The histogram in the left panel uses the scale of CG, the right one 
refers to the scale of ZW. The line at the top right indicates the median measurement error of 0.17dex. Overplotted as solid lines are MDFs, 
convolved by the respective errors. 
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B-V B-V 

Fig. 10. — Left panel: Distribution of metallicities across the RGB. Red squares depict the metal-poor RGB stars with [Fe/H]cG ^ —1.68, 
whereas blue circles refer to the metal-rich tail of the distribution ([Fe/H]cG > —1.68). Also indicated as solid lines are globular cluster 
fiducials from Sarajodini & Layden (1997) at the respective [Fe/H]cG, from left to right: M15 (-2.02), NGC6752 (-1.24), NGC 1851 (-1.03) 
and 47Tuc (—0.78). Representative color errorbars are indicated in the left panel. — Right panel: The same data points shown with Yonsei- 
Yale isochrones (solid lines; Yi, Kim, & Demarque 2003). From left to right, the isochrones have the following parameters: [Fe/H] = —2.3 
dex and 12.6 Gyr, —1.7 dex and 6.3 Gyr, —1.3 dex and 3.2 Gyr. The choice of parameters corresponds roughly to the three main episodes of 
an assumed star formation history of Carina. Comparing the location of the isochrones with the red giant branch fiducials on the left-hand 
side illustrates the effects of the age-metallicity degeneracy. 
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Fk;. 11. — CaT metallicities of the 437 red giants versus their clhptical radius in Carina (top). Nominal core and tidal radii are also 
denoted. The bottom left panel shows MDFs for three regions at different radii. These density distributions were convolved by the individual 
measurement errors in [Fe/H]. A weak trend of the MDF to become more metal-poor when proceeding outwards is visible and also reflected 
in the cumulative distribution (bottom right). 
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Elliptical Radius [arcmin] 

Fig. 12. — Cumulative number distributions for elliptical radii of the metal-poor population ([Fe/H]< —1.68, solid line) and the more 
metal-rich component (dashed line). The latter is more centrally concentrated. 




Fig. 13. — Fits of basic models ( modified simple closed box with outflows and prompt initial enrichment (PIE)) overplotted on the 
frequency distribution of the observed metallicities, where the dash-dotted line refers to the predicted distribution for Carina from Lanfranchi 
& Matteucci (2004, LM04) , scaled to match the total number of stars. The right panel is the analogous plot where the models were convolved 
by the observational errors. 



Table 1 

Globular Clusters used for calibration of the metallicity scale. 



Cluster 


a (J2000) 


6 {J2000) 


[Fe/HJzw^,ii97 


[Fe/HJcG,iJ97 


(m-M)y 


E(B-V) 




[Ca/FeJ 


References 


NGC 3201 


10 17 36.8 


-46 24 40 


-1.53+0.03 


-1.24 + 0.03 


14.21 


0.21 


14.74+0.07 


0.11 


1,2,3,4 


NGC4147 


12 10 06.2 


+ 18 32 31 


-1.77 + 0.04 


-1.50 + 0.06 


16.48 


0.02 


16.95 + 0.10 




1,2,3 


NGC 4590 


12 39 28.0 


-26 44 34 


-2.11 + 0.04 


-2.00 + 0.03 


15.19 


0.04 


15.75 + 0.05 


0.32 


1.2.3.4 


NGC 5904 


15 18 33.8 


+02 04 58 


-1.38 + 0.05 


-1.12 + 0.03 


14.46 


0.03 


15.13 + 0.05 


0.21 


1.2,3,4 



References. — 



(1) Harris (1996); (2) Rutledge et al. (1997); (3) Ferraro et al. (1999); (4) Carney (1996) 
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Table 2 
Observation log 



Date 


Field (Configuration) 


Total exp. time [sj 


21 Feb 2003 


Center 


1560 




NGC3201 


627 


22 Feb 2003 


Center 


9900 




NGC4147 


2400 


03 Mar 2003 


Center 


9497 


04 Mar 2003 


Center 


7800 




NGC 5904 


3600 


05 Mar 2003 


Center 


7800 




NGC 4590 


1800 




NGC 3201 


1200 


22 Dec 2003 


NE (Ic) 


14981 




NE (lb) 


12975 


23 Dec 2003 


Center (2b) 


13630 




Center (2a) 


13213 


24 Dec 2003 


SW (le) 


14401 




SW (Id) 


13172 


25 Dec 2003 


NW (le) 


14333 




NW (Id) 


13113 


26 Dec 2003 


NE (2b) 


14470 




NE (2a) 


13001 


27 Dec 2003 


SE (le) 


20488 




SE (Id) 


13333 


28 Dec 2003 


SW (2a) 


19469 




SW (2b) 


235 




NW (le) 


7933 


29 Dec 2003 


SW (2b) 


13421 




Center (le) 


13553 


30 Dec 2003 


Center (Id) 


13694 




Center (2a) 


13393 


21 Feb 2004 


SE (2a) 


19677 


22 Feb 2004 


SE (2b) 


15162 




NGC 4590 


4794 




NGC 3201 


3598 


23 Feb 2004 


NW (2a) 


10824 




NGC 4590 


6596 


24 Feb 2004 


NW (2b) 


10833 




NGC 3201 


7196 


25 Feb 2004 


NE (Id) 


23481 


26 Fob 2004 


NE (Ic) 


10884 


27 Feb 2004 


NE (2a) 


10931 


28 Feb 2004 


NE (2b) 


10635 



NOTE - Labels in brackets refer to different configurations 
on the fiber positioning plates, thus aiming at different targets. 
Only nights relevant for Carina observations are listed. 
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Table 3 

Observed fields in Carina and the calibration clusters 



Field 


a (J2000) 


6 (.J2000) 


Center 


06 41 36.8 


-50 57 58 


NW 


06 41 00.8 


-50 45 06 


NE 


06 42 51.1 


-50 52 50 


SW 


06 39 50.0 


-51 04 49 


SE 


06 42 05.3 


-51 10 00 


NGC 3201 


10 17 36.8 


-46 24 40 


NGC4147 


12 10 09.6 


+ 18 35 511 


NGC 4590 


12 39 28.0 


-26 44 33 


NGC 5904 


15 18 03.5 


+02 06 41 



Table 4 





Stars in common with Shetrone et al. 


(2003). 




Star ID 
(Mateo 1993) 


EIS-ID 


[Ca/FeJ [Fe/HJhighrcs 
(Shetrone et al. 2003) 


[Fe/HJcaT,ZW 

(This 


[Fe/HJcaT,CG 

work) 


M2 

M3 

M4 

MIO 

M12 


LG04a_003764 
LG04a_001419 
LG04a_001673 
LG04d.006644 
LG04a_007126 


+0.20 + 0.05 -1.60 + 0.02 
-0.10 + 0.06 -1.65 + 0.02 
+0.14 + 0.04 -1.59 + 0.02 
-0.02 + 0.05 -1.94 + 0.02 
+0.12 + 0.05 -1.40 + 0.02 


-1.67 + 0.09 
-2.00 + 0.08 
-1.65 + 0.09 
-2.16 + 0.08 
-1.33 + 0.11 


-1.43 + 0.09 
-1.84 + 0.08 
-1.41 + 0.09 
-2.04 + 0.08 
-1.01 + 0.11 



NOTE — Metallicities given for the Shetrone et al. (2003) stars are the weighted 
means from [Fel/H] and [Fell/H]. 
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Table 5 

sta.rs ileasi'eed l\ more than one obhermnc rtks. 





[rc/xij(7(;;!, nrst run 


re/ rl 


CGi second run 


L\ re/ n\cG 


[Fc/HJ(7(7, combined spectra 


ijLiU4a_uuuuoU 


1 /1 _i_ n 1 1 
— i.4y zt u.ii 




1 77 -U n 1 n 

1. { 1 it U.iU 


U.zo 


1 rcQ _i_ n 1 n 
— i.oo it U.iU 




1 /I c; _|_ n 1 /I 
— 1.40 zt U.i4 




1 /I c _i_ n no 
i.4o it u.uy 


u.uu 


1 /I r; _i_ n 1 n 
— 1.40 it U.iU 


ijLiU4a_UUU40l 


1 of; -L n 
— i.yo =t U.oD 




1 AK 1 n 1 
i.40 it U.iz 


— U.oU 


1 iri 1 n 1 Q 
— i.oi it u.io 


ijlorU4a_UUiU4Z 


1 tto 1 n 1 n 




1 Qfi 1 n 1 n 
i.oD it U.iU 


—U.i / 


1 ACi 1 n 1 n 
— i.4y it U.iU 


ijVjU4a_uUi i 1 1 


1 (^n 1 n 1^1 
— i.DU =t U.oi 




1 71^ 1 n 1 Q 
i. < it U.io 


U.io 


1 AO 1 n 1 Q 
— i.Dz it U.io 


ijLrU4a_UUi i4z 


ni J- n 1 




1 Ai 1 nil 
i.oi it U.ii 


— U.4U 


1 QK 1 nil 
— 1.00 it U. 1 1 


T r^fi/i '1 nni 1 vn 
ij^U4a_uui i / u 


1 AO 1 n 1 




i.Oo it u.io 


n no 
— U.uy 


1 AO -u n 1 n 
— i.Dz it U.iU 


T r^n/i nni oq/1 


n 01 -U n /I /I 
— u.yi =t U.44 




i.Oo it U.i / 


u.oz 


1 /I /I -U n 1 
— i.44 it u.iy 


J r^f\An nni ooq 
ijiorU4a_uuizyo 


n 70 J- n Qo 

— u. / y =t u.oz 




QO 1 n 1 
z.oz it U.iy 


L.OO 


no 1 n 1 Q 
— z.uy it U.io 


T r^n/1 nni qi q 


1 nA 1 n 01 
— l.UD =t U.zl 




OA 1 n 1 1 
z.z4 it U.14 


i.io 


1 Q J- n 1/1 
— z. lo it U. 14 


T r^n/io nni 
ijlafU4a_UUloD4 


1 Ai J- n 1 Q 
— i.Di =t U.io 




1 KT J- n no 

i.o / it u.uy 


— U.U4 


1 AC 1 n 1 n 
-i.Oo it U.iU 


J r^HAn nni Qon 
ijLiU4a_uuioyu 


1 AQ -U n 1 Q 

— i.Oo =t u.io 




1 fio -U n 01 
i.Dz it u.zi 


— U.UO 


1 QO -1- n 1 Q 

— i.oy it u.io 


T on/i nni a op^ 

ijLjU4a_UUi4zD 


1 00 -U n 
— i.yz zt u.o ( 




1 KQ -U n 1 1 
i.Oo it u.ii 


— U.o4 


1 AQ -U n 1 1 
— i.Do it U. 1 1 


T r^n/i '1 nr^i r^r^fi 


— 1.4U it u.iu 




1 /1 1 _i_ n no 
i.4i it u.uy 


U.Ui 


1 /1 1 _i_ n 1 n 

— 1.41 it U.IU 


T r^n/io nni i^c^q 


1 Q 1 1 1 r 

— i.oi it u.io 




1 in 1 n 1 n 
i.4U it U.IU 


— U.4i 


1 AO 1 n 1 Q 
— I.Dz it U. 1^5 


T r^n/i nni Ano 
ijLiU4a_UUlDUo 


1 ACi 1 n 1 1^ 
— i.4y it U.io 




1 A7 J- n no 
i.D/ it U.uy 


U.io 


1 cn 1 nil 
— i.oU it U.ii 


T r^n/io nni tqa 


1 en J- n 00 
— i.oU it yj.ZZ 




1 OQ 1 n 1 n 
i.yo it U.iU 


U.io 


1 7C 1 nil 
— i. /o it U.ii 


T r^n/i nni c7p; 


1 J- f\ ^ A 

— i.Oo it U.i4 




AQ -\- r\ ^ K 

z.4o it U.iO 


u.yo 


nn -U n 1 
— z.UU it u. iz 


T r^n/i nni coo 
ijL}rU4a_uuioyy 


1 _i_ n QO 
— i.Oo it u.oy 




ni -U n 1 
z.Ui it U.iz 


U.4D 


1 00 -u n 1 
— i.yz it u.iz 


T r^n/i nni 01 n 
ijLTU4a_uuiy iu 


1 KCi -U n 1 

— i.OU it U.iz 




1 Ai -U n 1 
i.oi it u.io 


nil 
u.ii 


1 f;a I n 1 .1 
— l.oD it U. 14 


T r^n/1 nni 017 

ijlarU4a_UUiyi / 


1 AK 1 n 1 Q 
— i.4o it u.io 




1 AA 1 n 1 Q 
i.oo it U.io 


n 01 
U.zi 


1 An 1 n 1 
— i.OU it U.iz 


T r^r\An nnoni^t=; 


1 QT 1 nil 
— L.of it u.ii 




1 _L n no 
i.oo it U.uy 


U.io 


1 Ko 1 n 1 n 
— i.oz it U.iU 


T r^n/i nnoi pc\ 
ijLTU4a_uUziDy 


1 Kr\ J- n nn 

— i.ou it u.uy 




1 1^0 J- n no 
l.oz it U.uy 


n no 
U.Uz 


1 1 j_ n 1 n 
— i.oi it U.IU 


T r^n/i nnoi ci 
ij^U4a_UUzioi 


1 /I /I _i_ n no 
— i.44 it u.uy 




1 Q7 -U n no 
i .0 ( it u.uy 


n n7 

— U.U / 


1 QO -U n 1 n 
— i.oy it U.IU 


ij^U4a_uuoyoz 


1 oc; -U n 1 
— i.zo it u.iy 




1 ao _i_ n 1 1 
i.Oo it U.ii 


n /I Q 
U.4o 


1 A/1 -U n 1 /I 
— 1.D4 it U. 14 


ij*orU4a_UU4i (i) 


1 on 1 n 1 Q 
— i.oU it U.io 




1 00 J- n no 
i.yz it U.uy 


n 1 
U.iz 


1 OQ 1 n 1 n 
— i.yo it U.iU 


T r^n/iv» nnoi^QQ 

ij(orU4D_UUZDC)0 


1 OK J- n OT 
— i.Zo it U.z / 




1 n J- n 11 
z.iU it U.ii 


U.oO 


1 C7 1 n 1 Q 
— i.o / it U.io 


T r^n/ii-v nni^on/1 
ljLrU4D_UUozy4 


1 00 J- n 01^ 
— i.zz it U.zo 




1 on 1 n OK 
l.yu it U.zo 


n i^Q 
U.Do 


1 c; 1 j_ n 1 
— i.oi it U.iy 


T r^r\Ar> nnnf^of; 
1j*^U4C_UUU0Z0 


1 01 1 n 1 n 
— i.yi it U.iU 




1 1 J- n no 
z.ii it U.Uo 


n on 
U.zU 


OA J- n no 
— z.Uo it U.uy 


T r^nyi/^ nnciT^n 

ijVjrU4C_UUU / iU 


1 yi/i _i_ n Ad 
— i.44 it u.4y 




1 OA -\- r\ ^ K 
i.z4 it u.io 


n on 
— U.ZU 


1 1 Q -U n 1 y1 
— i.iy it u.i4 


T r^n/ir. nnQnQc: 


1 1 -U n 1 Q 
— i. iz it u.io 




1 1 A -U n 1 1 
i.iO it u.ii 


n n/i 

U.U4 


n OA -U n 1 /f 
— U.yo it u. 14 




1 vn 1 n oi 
— i. / U it U.zi 




1 A1 j_ n on 
i.oi it U.zU 


n no 

—U.uy 


— I.OD it U.iz 


LG04c_004227 


-2.32 ±0.09 




2.53 ±0.10 


0.21 


-2.31 ±0.10 


LG04c_004308 


-2.08 lb 0.10 




2.13 ±0.08 


0.05 


-2.15 ±0.10 


LG04c_006477 


-1.25 lb 0.11 




1.75 ±0.09 


0.50 


-1.55 ±0.10 


LG04c_006479 


-1.51 lb 0.22 




1.75 ±0.09 


0.24 


-1.65 ±0.10 


LG04c_006573 


-1.52 lb 0.10 




1.42 ±0.09 


-0.10 


-1.45 ±0.10 


LG04c_006593 


-1.34 lb 0.20 




1.36 ±0.20 


0.02 


-1.39±0.11 


LG04c_006601 


-1.52 lb 0.31 




1.53 ±0.09 


0.01 


-1.54 ±0.10 


LG04c_006788 


-1.08 lb 1.07 




1.42 ±0.13 


0.34 


-1.38 ±0.13 


LG04c_007260 


-1.58 lb 0.09 




1.49 ±0.09 


-0.09 


-1.52±0.10 


LG04d_003625 


-2.01 lb 0.20 




1.60 ±0.34 


-0.41 


-1.68±0.18 


LG04d_004311 


-1.90 ±0.53 




2.14 ±0.28 


0.24 


-1.85±0.12 


LG04d_004565 


-1.65 ±0.69 




2.06 ±0.14 


0.41 


-2.05 ±0.14 


carl.tlOO 


-2.12 ±0.15 




1.90 ±0.08 


-0.22 


-2.01 ±0.11 


carl_tl02 


-1.62 ±0.10 




1.53 ±0.09 


-0.09 


-1.55 ±0.10 



Table 6 

Derived properties of member stars in Carina 



Star" 


a (J2000) 


5 (J2000) 


r['] 


V 


B-V 


EW[A] [Fe/H]fl97,w 


[Fe/H]ra7,cG 


cr [Fc/H] 


logp 


LG04a.00002 


06 


42 


59 


-51 


00 


06 


15.4 


19.50 


0.86 


3.81 


-1.82 


-1.62 


0.16 


-1.31 


LG04a.00008 


06 


42 


31 


-50 


59 


59 


10.6 


18.55 


1.13 


4.01 


-1.92 


-1.74 


0.09 


-1.31 


LG04a.00017 


06 


42 


26 


-50 


59 


43 


09.6 


19.71 


0.94 


2.72 


-2.11 


-1.99 


0.39 


-1.32 


LG04a.00030 


06 


41 


58 


-50 


59 


32 


04.8 


17.93 


1.27 


4.29 


-1.94 


-1.77 


0.10 


-1.27 


LG04a.00046 


06 


40 


52 


-50 


59 


12 


07.4 


19.56 


0.83 


2.30 


-2.27 


-2.18 


0.20 


-1.46 


LG04a.00064 


06 


42 


52 


-50 


58 


48 


13.4 


19.86 


0.84 


4.08 


-1.68 


-1.45 


0.33 


-1.31 


LG04a.00111 


06 


41 


04 


-50 


57 


41 


05.8 


19.57 


0.83 


2.68 


-2.15 


-2.03 


0.47 


-12.92 


LG04a.00115 


06 


42 


08 


-50 


57 


32 


05.2 


19.70 


0.90 


2.63 


-2.14 


-2.02 


0.11 


-1.30 


LG04a.00160 


06 


41 


56 


-50 


56 


51 


03.3 


19.29 


0.92 


3.17 


-2.05 


-1.91 


0.10 


-1.55 


LG04a.00189 


06 


41 


27 


-50 


56 


05 


03.6 


18.84 


1.02 


4.04 


-1.86 


-1.67 


0.21 


-1.29 



NOTE — This Table is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for 
guidance regarding its form and content, r and E W denote the elliptical radius and the CaT linestrength, respectively. [Fe/H] and 
the respective uncertainty a are given following eqs. 1,2. Finally, p denotes the targets' membership probability. 

The nomenclature is such th^t LG04a-d designates the ElS-fields targeting Carina (ElS-team, priv. comm.) followed by the 
number in the respective input catalog. 



